Examining the crossover from hadronic to partonic phase in QCD 
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It is argued that, due to the existence of two vacua — perturbative and physical — in QCD, the 
mechanism for the crossover from hadronic to partonic phase is hard to construct. The challenge 
is: how to realize the transition between the two vacua during the gradual crossover of the two 
phases. A possible solution of this problem is proposed and a mechanism for crossover, consistent 
with the principle of QCD, is constructed. The essence of this mechanism is the appearance and 
growing up of a kind of grape-shape perturbative vacuum inside the physical one. A dynamical 
percolation model based on a simple dynamics for the delocalization of partons is constructed to 
exhibit this mechanism. The crossover from hadronic matter to sQGP as well as the transition from 
sQGP to wQGP in the increasing of temperature is successfully described by using this model with 
a temperature dependent parameter. 
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The theory of strong interaction 



quantum chro- 



modynamics QCD has a complicated phase structure 
It has been shown by lattice gauge theory [2] that at low 
temperature and density there is a confined, chiral sym- 
metry breaking phase with hadrons as basic elements, 
while at high temperature and density is the deconfined, 
symmetry restoration phase with partons (quarks and 
gluons) as basic degree of freedom. It is found that at 
zero baryon density and high enough temperature, the 
transition from hadronic to partonic matter is of a cross- 
over type 3*3. 

Crossover is a gradual change of the system from one 
phase to the other without a definite transition point. 
The crossover between hadronic and partonic phases at 
high temperature has been firmly settled by lattice-QCD 
from thermodynamic argument, but really what happens 
in crossover; how does the system crossover from one 
phase to the other, are still open questions. 

In this respect, it is worthwhile to consider a simi- 
lar problem in QED — the Mott transition in atomic 
gas, where a neutral atom gas is transformed to electro- 
magnetic plasma through the ionization of atoms one by 
one. However, such a mechanism is inapplicable to QCD, 
because the crossover of a hadron matter to quark-gluon 
plasma through the decomposition of hadrons to quarks 
(anti-quarks) one by one contradicts color confinement. 
According to color confinement, isolated color object can- 
not exist in physical vacuum. The energy of an isolated 
color object in physical vacuum is infinite 

Unfortunately, most of the current models on the mar- 
ket, which claimed crossover have not taken this 
color confinement property of QCD into proper account. 
In these models the initial parton system evolutes and 
hadronizes gradually. The hadronization condition is 
either that a parton ceases to interact with other par- 
tons or that the total color interaction from a pair 
(or a three particle state) of quarks with the remaining 
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system vanishes [6|. In both cases the hadronization is 
carried out one by one. After a long time, the system in 
consideration is dominated by hadrons but there are still 
a few partons moving in the physical vacuum of hadronic 
matter Q, violating the confinement property of QCD. 

The difficulty of crossover in QCD lies in the fact that 
QCD has different vacua — perturbative and physical — 
with partons and hadrons as basic element, respectively. 
How to realize the transition from one vacuum to the 
other during the gradual change — crossover — between 
partonic and hadronic phases is a big challenge Q. In 
this letter we will discuss this problem and try to find a 
possible mechanism for the crossover in QCD compatible 
with the principle of color confinement. 

For this purpose let us turn to another kind of model — 
the geometrical bond percolation model fioj . This model 
has no dynamical prescription but is nevertheless enlight- 
ening in constructing a crossover model for QCD 

In the bond-percolation model, crossover is realized 
through the formation of clusters. In the lattice-version 
of the model, "a cluster is defined as a group of nearest- 
neighboring occupied sites that are linked by occupied 
bonds. The cluster size is the number of sites in a clus- 
ter' In our case, the sites are distributed continu- 
ously in the space and will be in the following referred to 
as cells. The cluster- formation could be understood as: 
originally, isolated cells are color singlet hadrons, when 
they are connected by bonds to form clusters, color can 
flow among them through bonds, and only the cluster as 
a whole keeps to be color singlet. 

Clusters could be of various sizes. The formation of 
an infinite cluster, which is defined in a finite system 
as a cluster extending from one boundary to the other, 
i.e. a cluster having, at least in one dimension, the sys- 
tem size, is taken as the appearance of a new constituent 
in the system. At that time there are still a lot of clus- 
ters of various-sizes, cf. Fig. 1(a), i.e. the system as a 
whole has not yet turned to the new phase. So we con- 
sider the appearance of an infinite cluster as the starting 
point of the crossover of the system to a new phase. The 
crossover process is completed when all the cells in the 



2 




FIG. 1: (Color online) Continuously distributed cells con- 
nected by bonds form clusters. In (a) the big cluster extend- 
ing from the left- to the right-boundary is an infinite cluster. 
In (b) all the cells are connected to an infinite cluster. 




Molecular- 
like 



grape shape perturbative vacuum [l5| of the nuclear size 
in 3-D, cf. Fig. 1(b). Meanwhile, there are still isolated 
bubbles of physical vacuum in the system [l^ . cf. the 
white regions in Fig. 1(b). The groups of colored parton 
are confined in the separate wells, or cells, inside the grape 
shape perturbative vacuum, being able to exchange color 
among themselves via quantum tunneling through the 
potential barriers. 

It could be expected that, the physical- vacuum bub- 
bles can not survive long in the environment of the grape- 
shape perturbative vacuum and will soon be transformed 
to perturbative too. The groups of colored parton can 
then move around in the whole system, resulting in a 
quark gluon matter possessing the property of perfect 
fluid, which is conventionally referred to as strongly cou- 
pled QGP, sQGP 

At still higher temperature the potential barriers be- 
tween neighboring wells will drop to zero and all the wells 
disappear. In this case all the parton-groups will be dis- 
assembled to independent partons and the sQGP turned 
to weakly coupled QGP, wQGP. 

Let us try to construct a simple model for realizing the 
above assumptions and arguments. Our aim is to answer 
the question: how to crossover in QCD. For simplicity, 
we take the initial system to be a nucleon gas. 

The model Hamiltonian for the 6 quarks in two near-by 
cells is assumed to be [l3 | 



FIG. 2: (Color online) Two possible forms of nucleon (or 
hadron) aggregation: gas-like and molecular-like. In both 
cases only the whole aggregation is color singlet. 



system are connected to a unique cluster, cf. Fig. 1(b). 
In this process no isolated color object will appear in 
physical vacuum, and thus no contradiction with QCD. 

Inspired by this achievement we propose our ba- 
sic assumption as the following. The aggregation of 
hadrons could be of two forms: gas-like or molecular-like, 
cf. Fig. 2 [l3[. In the first case all the hadrons in an ag- 
gregation are melted together, forming a big "bag" with 
quarks and gluons moving inside, which is usually re- 
ferred to as QGP droplet. In the second case, the hadrons 
(or cells) in a cluster, like the atoms in a molecular, in 
spite of being no longer color singlets, still keep their in- 
dividuality, i.e. still occupy separate space regions, con- 
nected with each other by bonds, forming a color singlet 
cluster. 

The cells, i.e. the separated space regions occupied by 
the colored parton-groups, are essentially "holes" dug 
by the groups of colored parton in physical vacuum [1]. 
These holes can be viewed as potential wells if we add 
energy as an extra ordinate. Partons can tunnel through 
the potential barriers between neighboring wells all over 
the color-singlet cluster, which is referred to as delocal- 



ization [14 1 



When the crossover process is completed, all the wells 
in the system are connected to a unique cluster, forming a 
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where Tcm is the center-of-mass kinetic energy. is 
the color interaction, which will be chosen as a square- 
confinement potential ~ —OcXi ■ Xji'^j when the 
quarks i,i belong to one and the same cell. When they 
belong to two nearby cells, the infinite potential in be- 
tween will drop down, forming a potential barrier, and a 



parametrization — 



— — af.A,; • A,- 



will be used. 



where /i is a model parameter. It turns out that as the in- 
creasing of pL the maximum distance 5o for delocalization 
increases, cf. Fig. 3 below. Since at higher temperature 
the quarks will be more free to move, i.e. 5*0 will be 
larger, we require pL to be an ascending function of T. 
From dimensional consideration we assume /i oc T^. 

In doing variational calculation the trial wave function 
of the two-cell system in adiabatic approximation is cho- 
sen to be an antisymmetric six-quark product state [l3 | 



= 1 1=4 J 00 

(2) 

where A is the anti-symmetrization operator, which per- 
mutes quarks between the two cells; [• • • ]oo means that 
the spin, isospin and color of the two cells are cou- 
pled to a particular color singlet state with total spin 
and isospin equal zero. For the orbital motion, we 
have the left (right) single-quark orbital wave function 
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where are the ceh-centers, S is the distance between 
the two cells. 6 is a baryon-size parameter. Delocal- 
ized orbit is defined as = "^['/'l + i^nif) — 

jf [£4'l+4'r] : where e is a variational parameter character- 
izing the degree of delocalization and N a normalization 
factor. At each separation S, e is determined by min- 
imizing the energy £(8) = ^'^^^(j.'^ffi?" ■ The model 
parameters are: TOu = Wd — 313 MeV, b = 0.603 fm, 
flc = 101.14 MeV/fm^, while fi is leaving as a free pa- 
rameter. 

It turns out that the model gives a reasonable amount 
of delocalization. When the two cells are close together, 
i.e. S is small, the delocalization is large e=l, but when 
the two cells separate and S increases to a certain dis- 
tance 5*0, the delocalization degree e suddenly drops to 
zero. The dependence of 5*0 on parameter fi is shown in 
Fig. [3 
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FIG. 3: The fi dependence of So determined by the dynamics 
of our model. Small pad shows a local part which is useful in 
determining the critical values of parameter ^. 

5*0 provides us a distance within which the quarks have 
equal probability to be simultaneously in the two cells, 
i.e. within which bond is likely to be formed. We will 
use it to construct a bond-percolation model to realize 
the crossover in QCD. For simplicity, we will consider 
the system to be a 2 dimensional static one. 

The initial configuration consists of 2 x 197 cells, which 
are small spheres (circles in 2-D) of hard-core radius re = 
0.1 fm distributed randomly in a big sphere (circle in 2- 
D) of radius i? = 7 fm. A cell with center departing from 
the center of the big sphere (circle) farther than R — r^ , 
is considered as a boundary cell. 

The percolation procedure is as follows: 

1) Randomly select a cell a as a mother cell. 

2) Define the cells with |r — rQ,| < 5*0 as bond- candidate 
cells, randomly select 3 of them to form bonds connected 
to the mother cell a, and call the latter as daughters. 



If the number of candidate-cells is less than 3, then the 
number of daughters is equal to the candidate number. 

3) For every daughter of cell a randomly select 2 
bond-candidate cells of it to form bonds, and call them 
granddaughters. Repeat the procedure until no bond- 
candidate cell can be found anymore. 

4) Then choose another cell (3 from the left unbounded 
cells as another mother cell, and repeat the procedure 
starting from step 2. 

5) Repeat step 4 until all the cells are exhausted. 

In this way, every cell is assigned to a cluster. In each 
cluster, find the boundary cells if any, and calculate the 
distance between every two boundary cells. Denote the 
maximum distance by d. A cluster with d > y/2R is 
taken as an infinite cluster. 
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FIG. 4: (Color online) Numerical results Poo (full circles) 
and TVs/lOO (open triangles) versus So for 10 000 events. The 
dashed lines correspond to the 2 thresholds, i.e. Tc the place 
where Poo starts to increase from zero, T^' where the mean 
number Na of cells outside an infinite cluster tends to zero. 

A configuration with all the cells assigned to clusters 
is called an event. Generate N events, and determine 
in each event whether there are infinite cluster (s). Sup- 
pose M is the number of events with infinite cluster(s), 
the probability for the appearance of event with infinite- 
cluster(s) can be expressed as Poo — limAr^oo M/N . The 
dependence of Poo on 5*0 is plotted in Fig. 4 as full cir- 
cles. We see that at a certain value 5*0 = Sc, Poo starts 
to increase from zero, and the system starts to crossover 
to a new phase. The corresponding value of is = ^c- 

Let Ns be the number of cells outside of an infinite 
cluster in an event. The mean versus S'o is plotted in 
Fig. 4 as open triangles. The point where TVs = marks 
the accomplishment of crossover and will be denoted by 
S'o = Sc'. The corresponding value of /i is /ic', cf. the 
dashed lines in Fig. 4 and in the small pad of Fig. 3. 

The starting and ending points Hc and ^c' of crossover 
can be extracted from Fig's. 3, 4. Assuming oc we 
get the corresponding temperatures: Tc' ~ 1.39 Tc. 
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It is interesting to see from Fig. 3 that as fi increases 
further from fi,.', at a certain point fic", Sq increases 
sharply to infinity. An infinite 5*0 means that partons 
can flow from one well to the other no matter how far 
away they are. This means that at this point the poten- 
tial barriers have dropped to zero, and there is no well in 
the vacuum anymore. Partons can move freely inside the 
whole system, and the sQGP turns to wQGP. Extract- 
ing the point /j,c" for this transition from the figure and 
assuming again esc T^, we get Tc" ~ 1.98Tc p^ . 

In this letter a possible mechanism for the crossover 
in QCD, compatible with the principle of color confine- 
ment, is proposed. The way for vacuum-change during 
the gradual crossover between the two phases is via the 
following steps. 

1) As the increasing of temperature, partons in hadrons 
start to be delocalized, tunneling between neighboring 
hadrons, and change the latter to colored parton- groups 
located in potential wells. 

2) The wells connected by tunnels form color singlet 
clusters, which grow up in the physical vacuum, and even- 
tually become a big cluster of the system size. Inside this 
cluster is a grape-shape perturbative vacuum, while out- 
side of it is physical vacuum bubbles, which will change 
to perturbative as the further increase of temperature, 
and the colored parton-groups in the wells can then move 
around, forming a fluid-like matter — sQGP. 

3) The potential barriers will drop to zero at a higher 
temperature and all the wells disappear, the sQGP is 



then turned to wQGP. 

Last but not least, a dynamical percolation model is 
constructed for the first time, which is superior to the 
purely geometrical percolation models. 1) Our model 
has dynamical foundation and the cells and bonds have 
definite physical meaning. 2) We use a maximum length 
Sq for bond-formation, obtained from dynamical calcula- 
tion, and construct 3 bonds randomly within this length. 
This is to substitute the bond-formation probability p 
put in the geometrical models by hand. 3) Our model 
has a temperature-dependent parameter fi and is able 
to deal with the crossover as well as the transition from 
sQGP to wQGP in the increasing of temperature, while 
the geometrical percolation models are able to treat only 
the effect of the increasing of density. 

In the proposed scenario the crossover is accomplished 
within a temperature range from Tc to Tc' . It could be 
expected that at some energies lower than 200 GeV, the 
produced dense matter may be at a mediate stage of 
crossover, which in our scenario consists of an infinite 
cluster accompanied by a large number of other clusters, 
cf. Fig. 1(a). To study the properties of such a structure 
of matter theoretically and experimentally is worthwhile, 
which will deepen our understanding on QCD phase di- 
agram and critical dynamics. 
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